AO  A 0 7 2 7 1 4 


AFAL-'TR-79-1115,  Vol  I 


OPTICAL  FRINGE  ANALYSIS  TECHNIQUES 


Dr.  T.  L.  Williamson 
C.  R.  Lane 


Electro-Optica  Technology  Branch 
Electronic  Technology  Division 


r>  TECHNICAL  REPORT  AFAL-TR-79-1115,  Vol  I 


tLipinai  Report  for  Period  1 July  1975  - 30  September  1978 


Approved  for  public  release,  distribution  unlimited 


AIR  FORCE  AVIONICS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


NOTICE 


When  Government  drawings,  specifications , or  other  data  are  used  for  any  pur- 
pose other  them  in  connection  with  a definitely  related  Government  procurement 
■operation,  the  United  States  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever ; and  the  fact  that  the  government  way  have  formulated , 
furnished , or  in  any  way  supplied  the  said  drawings,  specifications , or  other 
data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner  licen- 
sing the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any 
way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS) . At  NTIS,  it  will  be  avail 
able  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


i2  rUL** 


CHARLES  RICHARDJAHE,  Engineer 

LiecCPo-OptlcaTechniques  and 
Applications  Group 
Electro-Optics  Technology  Branch 


FOR  THE  COMMANDER 


KENNETH  R.  HUTCHINSON,  Chief 
Electro-Optical  Techniques  and 
Applications  Group 
Electro-Optics  Technology  Branch 


RONALD  F.  PAULSON,  Acting  Chief 
Electro-Optics  Technology  Branch 
Electronic  Technology  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization  please  notify 
AF AL/DHO-2  , W-PAFB , OH  45433  to  help  us  maintain  a current  mailing  list”. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  se- 
curity considerations , contractual  obligations,  or  notice  on  a specific  document 


UNLLAt>blf  IhU 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wt »n  Of  Enltttd) 


READ  INSTRUCTIONS 


REPORT  DOCUMENTATION  PAGE 


BEFORE  COMPLETING  FORM 
J RECIPIENT'S  CATALOG  NUMBER 


GOVT  ACCESSION  NO, 


AF AL-TR- 79-1115,  Vol  I-- 


Final  Report 


OPTICAL  FRINGE  ANALYSIS  TECHNIQUES 


S CONTRACT  OR  GRANT  NUMBERS •) 


r — 

T . L. /Williamson 


10  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A UfORK  UNIT  NUMBERS 


» PERFORMING  ORGANIZATION  N AUl  ANU  ABB 

AIR  FORCE  AVIONICS  LABORATORY 

AFAL/DHO-2 

WPAFB,  OHIO  45433 


62204F 


12  RCPOUT 

Jun#  f§79 


n.  CONTROLLING  OFFICE  NAME  AND  AOORESS 


15-  SECURITY  CLASS,  (el  thla  report) 


MONITORING  AGENCY  name  A ADDRESS*”//  dll  latent  from  Controlling  Otlice) 


UNCLASSIFIED 


15a.  OECLASSI  FI  CATION /DOWNGRADING 
SCHEDULE 


U DISTRIBUTION  STATEMENT  (ol  thla  Report) 


Approved  for  public  release;  distribution  unlimited 


17  DISTRIBUTION  STATEMENT  (ot  (hrn  mbatrect  entered  In  Block  20.  if  ditleienl  from  Roporl) 


»•  supplementary  notes 


*t  KEY  WOROS  (Continue  on  reeerae  aid*  II  neceaaary  end  Identity  by  block  number  j 


Laser 

Interferometry 

Specklcgram 


ABSTRACT  (Continue  an  reveree  aide  II  neceeemry  end  Identity  by  block  number ) 


The  use  of  specklegrams  for  performing  certain  types  of  NDI  and  NDT  has  been 
known  for  some  time  yet  a method  for  extracting  the  data  quickly  and  cheaply 
has  eluded  specklegram  users.  This  report  discusses  a new  technique  for  per- 
forming the  data  reduction  of  the  large  amounts  of  data  obtained  through 
speckle  inter ferogr ams . This  new  technique  correlates  an  optical  mask  on  a 
fringe  pattern  to  generate  a single  line  scale  proportional  to  the  frequency 
of  the  cosine  fringe  pattern  when  Integrated  orthogonally  to  the  fringe  line. 


UNCLASSIFIED 

SECURITY  CLARIFICATION  OF  THIS  PAOE 


!•  Entered) 


TABLE  OF  CONTENTS 


SECTION  PAGE 

I INTRODUCTION  1 

II  SPECKLEGRAM  3 

III  SPECKLEGRAM  FRINGE  ANALYSIS  9 

IV  EXPERIMENTAL  VERIFICATION  16 

V CONCLUSIONS  19 

1 


SECTION  I 


INTRODUCTION 

Considerable  interest  has  been  kindled  over  the 
past  several  years  In  applying  holographic  Interferometry 
and  specklegram  Interferometry  to  the  areas  of  NDT  (non- 
destructive testing)  and  NDI  (non-destructive  Inspection). 
The  results  of  this  application  are  usually  less  success- 
ful than  desired.  The  primary  reason  for  this  Is  that  the 
interferometry  is  so  sensitive  that  the  analyst  is  over- 
whelmed with  data  and  must  spend  considerable  time  going 
over  this  data  to  deduce  what  It  all  means.  It  Is  .thus, 
essential  that  techniques  be  developed  to  reduce  the  fringe 
patterns  to  a lesser  quantity  of  more  meaningful  data. 

The  potential  for  automated  reduction  of  the  speckle- 
gram  fringe  data  has  been  demonstrated  and  Is  presently 
being  pursued  under  a joint  AFFDL/AFAL  program.  This 
program  utilizes  conventional  digital  processing  techniques 
to  perform  the  data  reduction.  The  required  system  for 
this  approach  Is  very  complex  and  expensive.  A technique 
Is  needed  that  Is  less  complex  and  hence  requires  less 
expensive  equipment  for  Implementation.  The  following 
discussion  Is  the  result  of  a first  attempt  to  take  a more 
general  look  at  the  problem  with  the  goal  of  a more  specific 
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and  simple  solution.  This  first  look  has  resulted  In  some 
promising  approaches  to  the  solution  of  the  specklegram 
fringe  analysis  using  optical  techniques  to  do  the  majority 
of  the  work.  This  results  In  lessening  the  requirements 
for  data  storage  and  hence  less  expensive  computer  equip- 
ment. 


SECTION  II 


SPECKLEGRAM 

The  soecklegram  under  consideration  Is  a photographic 
recording  of  a double  exposure  of  an  object.  The  speckle 
Is  the  result  of  using  laser  Illumination  on  an  object 
with  a diffuse  surface.  If  the  object  has  undergone  a 
surface  displacement  between  the  exposures  either  due  to 
external  forces  or  gross  displacement,  this  Information 
will  be  recorded  In  the  soecklegram.  The  ensuing 
discussion  Involves  the  recording  and  subsequent  readout 
of  this  displacement  Information. 

The  speckle  Is  In  effect  a random  sampling  function 
whose  sample  size  Is  small  compared  to  the  areas  under 
Investigation.  The  size  of  the  soeckle  can  be  controlled 
by  the  recording  optics.  The  speckled  Image  of  an  object 
can  be  represented  by 

lUj.yj)  Kx-xj.y-yj)  0) 

x1.Yj 

where  I(x,y)  represents  the  unspeckled  Image  of  the  object 
and  the  delta  function  represents  the  Individual  speckle. 
The  delta  function  Is  used  since  the  speckle  Is  small 
compared  to  the  minimum  resolution  of  the  Image  that  Is 
being  Investigated. 
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I(xk.y] ) S(x-xk,y-y1 ) 


(2) 


where  the  summation  Is  now  taken  over  the  second  set  of 
speckles.  For  small  surface  displacements,  the  positions 
of  the  speckles  on  the  second  Image  are  related  to  those 
of  the  first  by 


\ ■ xi  * 1 yi  ■ vj  * 


(3) 


where  Ax  and  Ay  represent  the  displacement  In  the  two 
coordinate  directions.  The  field  for  the  second  exposure 
can  then  be  written  as 

^ AX.Yj+  ay)  Mx  -x^  Ax  ,y  -y^-  Ay)  (4) 

The  photographic  recording  process  Is  a square  law  process 
resulting  In  a total  exposure  for  a unit  time  given  by 

E(x,y)  * 3 x1  ,yj  ) Mx  -x1 , y -yj  ) j + 

X1’yj 

jj(xi+Ax,yJ+  Ay)  ^(x-xf-  Ax.y-yj-  Ay)J  J(5) 
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The  summation  has  been  taken  outside  of  the  brackets 
since  the  cross  product  of  delta  function  Is  zero  unless 
the  arguments  are  equal . 

Since  the  specklegram  Is  being  used  only  for  a point 
to  point  comparison  of  the  surface  displacement,  the 
area  of  Interest  at  any  one  time  Is  small  compared  to  the 
detail  of  the  Image.  For  this  reason,  the  Image  function 
can  be  considered  a constant  over  each  area  of  Interrogat- 
ion. Thus,  over  a small  area  the  exposure  function  can 
be  represented  by 


E(x,y)  « K ( 2»  ( x-xi  ,y-yj))2+  (S(x-xi -6x  ,y-yj-&y) ) 

X1’yj 


(5) 


where  K Is  the  constant  value  for  the  Image  function  over 
the  area.  The  exposed  film  consists  of  a collection  of 
delta  function  type  sampling  functions  whose  positions  are 
determined  by  the  relationship  of  equation  6.  The 
granular  nature  of  the  photographic  emulsion  dictates 
that  either  the  silver  particle  exists  at  a given  location 
or  It  doesn't,  therefore,  the  squaring  of  the  delta 
functions  has  no  significance  In  the  resultant  Image  that 
Is  recorded.  The  resultant  recorded  Image  In  the  small  area 
can  be  represented  by 


E(x,y)  • K y 


( & ( x - x j ,y  ~y  j ) ♦ SCx-x^Ax.y-yj-Ay))  (7) 
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The  resultant  specklegram  will  then  consist  of  an 
array  of  speckle  sampling  functions  (represented  by  the 
delta  functions)  whose  positions  reflect  the  displacement 
that  occured  between  exposures.  Each  of  these  speckles 
will  be  recorded  In  the  form  of  an  absorptive  silver 
particle  or  collection  of  silver  particles.  It  will 
be  assumed  for  simplicity  that  the  film  will  be  completely 
absorptive  at  a speckle  position.  This  assumption  does 
not  affect  the  form  of  the  result,  but  rather  only  the 
multiplicative  constant.  The  transmission  of  the  film 
can  then  be  written  in  the  form 


T(x,y)  = 1 - ( Mx-x.  »y-yj  ) + 5(x-xi -&x  .y-yj-^y) ) 


Xi*yj 


3 1 -(S(x,y)  + Mx-fcx, y-Ay))  * Mx-Xj  .y-y^  ) 


«y. 


(8) 


where  * represents  the  convolution  operation. 

The  specklegram  is  analyzed  by  illuminating  each 
Individual  area  by  a small  collimated  laser  beam  and 
observing  the  resultant  far  field  pattern.  The  amplitude 
of  the  field  existing  in  the  far  field  region  of  the  area 
being  analyzed  will  be  the  Fourier  transform  of  the  field 
existing  at  the  exit  aperture  of  the  film.  If  the  field 
used  for  the  Illumination  is  of  unit  amplitude  then  the 
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exit  field  will  be  represented  by  the  expression  of 
equation  8.  The  field  In  the  far  field  region  will  then 
be  obtained  by  taking  the  Fourier  transform  of  this 
expression. 

Eff(x1 ,y1 ) = ^ [ T( x ,y ) ] (9) 


On  performing  the  Fourier  transform  term  by  term  the 
field  that  exists  In  the  far  field  can  be  represented 
by  the  expression 

Eff  (Xl  ,yi ) * 5(x1>yi)  + 2e-l(Axxl+  > 


Cos[(AxXj+  ^yy1 )/2] 


S(x1 - x ,j  ,y, -yj) 


(10) 


The  first  term  represents  the  "D-C"  or  transmitted  beam 
and  the  second  term  Is  a set  of  speckled  cosine  fringes. 
The  period  of  the  cosine  fringe  pattern  represents  the 
displacement  that  occured  between  exposures  and  Is  the 
desired  parameter  to  be  measured.  The  series  of  delta 
functions  Is  the  transform  of  the  series  that  existed  In 
the  exit  aperture  and  Is  merely  the  speckle  that  appears 
In  the  observed  pattern.  The  phase  term  will  not  be 
of  Importance  for  the  cases  when  the  Intensity  of  the 
far  field  Is  observed  rather  than  It's  amplitude.  A 
typical  far  field  Intensity  pattern  Is  shown  In  Figure  1. 
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Several  observations  can  be  made  at  this  point  which 
could  prove  to  be  of  great  benefit  in  establishing  any 
generalized  techniques  of  specklegram  analysis.  The 
speckle  will  be  present  at  all  stages  of  the  analysis 
unless  averaging  techniques  are  applied  to  eliminate  it. 
It's  size  can  be  controlled  through  proper  choice  of  both 
the  construction  and  analysis  optics  so  that  it's  effect 
may  be  minimized  in  the  analysis.  The  far  field  pattern 
is  not  general,  but  rather  a series  of  cosine  fringes. 

For  this  reason,  a generalized  analysis  of  the  pattern  is 
not  required.  The  desired  information  can  be  obtained 
from  only  the  fringe  spacings  and  their  orientation. 


FIGURE  1 : Typical  far  field  intensity  pattern. 


SECTION  III 


SPECKLEGRAM  FRINGE  ANALYSIS 


In  this  section  a technique  of  fringe  analysis  will 
be  discussed  that  relies  on  the  knowledge  of  the  form  of 
the  far  field  pattern  that  was  obtained  In  the  preceedlng 
section . 

Since  It  was  determined  In  the  preceedlng  section 
that  the  fringe  system  resulting  from  a specklegram  Is, 
or  at  least  can  be  made  to  be,  a series  of  speckled 
cosine  fringes,  the  following  analysis  applies.  Since 
the  far  field  pattern  Is  an  even  function,  the  Fourler- 
coslne  transform  can  be  used  to  determine  the  relationship 
between  fringe  pattern  and  the  speckled  double  Image.  In 
the  same  manner  the  Fourier  cosine  transform  can  be  used 
to  analyze  the  fringe  pattern  Itself.  The  fringe  pattern 
is  of  the  form  A Cosk0x  , where  A Is  the  arbitrary  ampl- 
itude and  kQ  Is  the  radian  frequency  of  the  pattern.  Ifi 
fact  If  the  expression  for  the  total  fringe  distribution 
Is  expanded  In  a cosine  series  only  one  term  would  have 
a non-zero  coefficient,  that  corresponding  to  the  proper 
spatial  frequency.  The  general  form  for  the  coefficients 
of  the  cosine  series  expansion  would  be 
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a*  " / 


ACos(kQx)Cos(kxx)  dx 


(11) 
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This  operation  can  easily  be  Implemented,  since  It  consists 
merely  of  multiplying  the  fringe  pattern  by  a cosine 
function  and  Integrating  over  the  variable  x . From 
this  It  Is  evident  that  If  a mask  were  constructed  with 
transmission  given  by  Cos(kj(y+b))  then  the  product  of 
the  incident  fringe  pattern  (when  properly  aligned)  and 
this  transmission  function  would  be  the  integrand  of 
equation  11.  The  integral  of  the  expression  over  the  x 
variation  would  be  a function  of  y and  would  represent  the 
coefficient  of  the  correspondi ng  cosine  term  in  the 
expansion.  The  expression  would  be 


ax(y)  = 


/ 


ACos { kQx )Cos ( k^ (y+b ) ) dx 


(12) 


This  expression  would  have  a non-zero  value  only  when 


k^y+b) 


(13) 


Thus,  the  y value  for  which  a non-zero  value  was  obtained 
would  be  a measure  of  the  fringe  spacing. 

In  reality  the  aperture  over  which  the  fringe  pattern 
exists  Is  not  Infinite,  resulting  in  the  expression 

,r> 

ACos(kQx)Cos(kj (y+b) ) dx  (14) 


a„(y> 


f 

n» 

/ 


P0(x)ACos(kox)Cos(k1 (y+b))dx 
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where  2D  Is  the  aperture  diameter  In  the  x direction  and 
PD(x)  takes  on  the  value  of  1 for  -D  < x<  D and  Is  zero 
elsewhere.  The  effect  of  this  non-infinite  aperture  Is 
that  Instead  of  a delta  function  occurlng  at  the  value  of 
y for  which  equation  13  holds  the  function  a x ( y ) will  be 
the  Fourier  transform  of  PD(x)  and  will  be  centered  at 
the  previous  value  of  y.  This  extends  In  the  y direction 
since  the  y variation  Is  a measure  of  the  x spatial 
frequency. 

ax(y)  - 2 S1n(k^D(y»b) ) (15) 

M ( v+b) 

This  function  will  determine  the  ultimate  resolution  of 
the  technique  as  a function  of  the  aperture  size  or  the 
number  of  fringes  analyzed. 

In  reality  the  transmission  of  the  mask  does  not 
have  to  be  a pure  cosine  function.  It  Is  only  required 
that  It  have  a fundamental  cosine  frequency  sufficiently 
strong  that  the  location  of  It's  maximum  can  easily  be 
determined  on  the  y-axls.  Therefore,  even  a square  wave 
variation  would  be  adequate.  The  spatial  variation  of 
such  a transmission  mask  would  appear  as  shown  In  Figure 
2. 

The  product  of  the  mask  funct  on  of  figure  2 with 
fringe  system  of  equally  spaced  lines  Is  shown  In  Figure 
3.  The  resultant  Moire'  fringe  systems  gives  a visual 


measure  of  the  correlation  that  Is  talcing  place.  The 
pattern  In  this  form  would  not  be  easy  to  measure  with 
an  E-0  device.  Figure  4 shows  the  function  that  Is 
obtained  when  the  function  In  figure  3 Is  Integrated 
over  the  x coordinate  (horizontal).  This  function  not 
only  readily  shows  the  location  of  the  y value  for  maximum 
correlation  but  also  shows  the  shape  of  the  aperture 
Induced  spread  function.  The  fringe  patterns  as  well 
as  the  Integration  In  these  figures  were  simulated  on  the 
Hewlett-Packard  9820  calculator. 

The  preceedlng  discussion  has  outlined  a rather 
simple  approach  to  determining  the  spatial  frequency  of 
the  cosine  fringe  pattern . The  procedure  used  does  not 
require  the  Implementation  of  a general  Fourier  transform 
operation.  As  will  be  shown  In  the  following  section 
this  technique  Is  easily  Implemented  optically. 


SECTION  IV 


EXPERIMENTAL  VERIFICATION 

In  order  to  demonstrate  the  concept  of  fringe 
analysis  discussed  In  the  previous  section  optically, 
a reference  specklegram  was  reconstructed  giving  the 
speckled  cosine  fringe  pattern  shown  In  Figure  5.  This 
fringe  system  was  Imaged  onto  a standard  mask  (Edmund 
Scientific  Moire'  Pattern  #4)  having  a square  wave 
transmission  characterl stl c whose  spacing  varied  propor- 
tionally to  the  radius.  The  resultant  product  pattern 
Is  shown  In  Figure  6.  Note  that  the  Moire'  fringes  can 
be  seen  In  the  same  form  as  those  of  Section  III. 


* * • _ \»  » . ^ 

ii; 


FIGURE  5 : Cosine  Fringe  System 


FIGURE  6 : Product  of  Cosine  Fringe  System  with 
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Transmission  Mask. 

The  resultant  field  Is  then  integrated  along  the 
direction  of  the  fringe  Dropoqation  by  means  of  a 
cylindrical  lens.  The  Intensity  pattern  that  is  obtained 
will  then  give  a measure  of  the  position  in  the  other 
coordinate  direction  where  the  correlation  between  the 
mask  and  the  fringe  oattern  Is  a maximum.  This  position 
can  then  be  related  to  the  spatial  frequency  of  the  mask 
at  that  position  and  the  resultant  fringe  soaclng  obtained. 
Figure  7 shows  the  result  of  this  Integration  on  the 
nattern  of  Figure  6. 
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FIGURE  7:  Resultant  Correlation  Function  Relating 
the  Snatial  Frequency  of  the  Fringes  to 
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Those  of  the  Mask. 

The  correlation  of  Figure  7 was  obtained  by  taking 
the  intensity  Drofile  of  the  image  on  an  image  analyzer. 

This  experiment  has  demonstrated  that  the  fringe 
analysis  can  be  implemented  simnlv  through  the  use  of 
a transmission  mask  and  a cylindrical  lens. 
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CONCLUSIONS 

The  oreceedlng  discussion  has  demonstrated  that  there 
are  techniques  of  analyzing  the  fringes  obtained  from 
specklegrams  that  are  easily  Implemented  optically. 

It  Is  not  claimed  that  the  technique  discussed  In  the 
preceedlng  sections  Is  the  only  or  the  best  solution  to 
the  problem,  but  rather  that  It  Is  one  solution  and  It 
does  work. 

The  procedure  that  has  been  outlined  takes  advantage 
of  the  known  characteristics  of  the  fringe  pattern  that 
Is  being  Investigated.  The  important  features  are  that 
the  fringe  pattern  Is  a cosine  function  and  that  It  Is 
always  centered  on  the  D-C  spot,  or  at  least  can  be  made 
to  be  centered.  Tnese  features  permit  the  analysis  to  be 
made  through  a rather  simple  correlation  technique  with 
one  transmission  mask.  The  resultant  correlation  function 
could  then  be  measured  by  a detector  array  to  give  a 
completely  automated  readout.  The  fringe  pattern  from 
the  specklegram  could  be  rotated  by  an  optical  system  using 
a dove  prism  before  correlation  to  obtain  the  angular 
Information . 


Before  the  techniaues  discussed  here  could  be 


implemented,  it  will  be  necessary  to  do  a more  thorough 
analysis  of  the  procedures  to  optimize  the  techniques. 

Some  of  the  areas  that  need  further  consideration  are; 

1.  Perform  an  error  analysis  of  the  procedure  to 
determine  the  accuracy  as  a function  of  the  number  of 
fringes  available  in  the  aperture  and  the  fringe  contrast. 

2.  Optimize  the  specklegram  recording  procedure. 

3.  Develop  techniques  of  measuring  the  fringe  spacing 
for  maximum  correlation  and  it's  associated  angle. 

4.  Determine  the  potential  of  combining  the  correlation 
and  integration  steps  into  a single  holographic  element. 
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